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^ ■ Abstract 
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. couplings, in nonrelativistic lattice QCD are determined to leading order in 



5 



X 



and radiative corrections to two important couplings, the so-called kinetic 



tadpole-improved perturbation theory. The scales at which to evaluate the 
running QCD coupling for these quantities, except the wave function renor- 



Oh! 

. malization, are obtained using the Lepage-Mackenzie prescription. When the 



bare quark mass is greater than the inverse lattice spacing, the kinetic cou- 
pling corrections are roughly 10% of the tree-level coupling strengths; these 
corrections grow quickly as the bare quark mass becomes small. A need for 
computing the two-loop corrections to the energy shift and mass renormal- 
ization is demonstrated. 
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I. INTRODUCTION 



Nonrelativistic lattice QCD (NRQCD) is an effective field theory formulated to reproduce 
the action of continuum QCD at low energies The NRQCD Lagrangian includes 

interactions which systematically correct for relativity and finite-lattice-spacing errors. To 
complete the formulation of NRQCD, the coupling strengths of its interactions must be 
determined. One possibility is to treat the couplings as adjustable parameters and tune them 
to fit certain experimental data; however, this tuning significantly reduces the predictive 
power of NRQCD simulations and is very costly and difficult. A better alternative is to 
compute the couplings in terms of the fundamental QCD coupling as = g^/ (47r) and the 
bare heavy-quark mass M using perturbation theory. This is done by evaluating various 
scattering amplitudes in both QCD and lattice NRQCD and adjusting the couplings until 
these amplitudes agree at low energies. Since the role of these couplings is to compensate for 
neglected low-energy effects from highly-ultraviolet QCD processes, one expects that they 
may be computed to a good approximation using perturbation theory, provided that ciAqcd 
is small, where a is the lattice spacing and Aqcd is the QCD scale parameter. 

In this paper, the heavy-quark propagator in NRQCD is calculated to leading order 
in tadpole- improved perturbation theory and the on-mass-shell dispersion relation for the 
heavy quark is obtained. The mass renormalization, a shift in the zero point of energy, 
and the radiative corrections to the couplings of two important NRQCD interactions are 
determined by matching this dispersion relation to that of continuum QCD to 0{v^), where 
V is the expectation value of the heavy-quark velocity in a typical heavy-quark hadron. 
The wave function renormalization is extracted from the on-shell residue of the perturba- 
tive propagator. These parameters are needed for high-precision numerical simulations of 
quarkonium and heavy-light mesons. Results are obtained using the standard Wilson action 
and a tree-level 0(a^)-improved action for the gluons, and using two versions of the NRQCD 
action: one containing all 0(f ^) spin- independent and all 0{v^) spin- dependent corrections 
to the leading kinetic term with tree-level removal of cutoff errors, the other including only 
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0(f ^) corrections with no tree- level removal of lattice-spacing errors. 

In addition, the scales at which to evaluate the running coupling for these quantities, 
except the wave function renormalization, are obtained for the first time using the Lepage- 
Mackenzie prescription. Choosing a value for aAqco from recent simulation measurements 
of the QCD coupling, numerical estimates of the energy shift, mass renormalization, and 
the radiative corrections to the two so-called kinetic couplings are also obtained. When 
aM > 1, where M is the bare quark mass, the radiative corrections to the kinetic couplings 
are roughly 10% of the tree-level coupling strengths; these corrections grow quickly as aM 
becomes small. A range of estimates for the energy shift and mass renormalization are 
obtained after circumventing certain defects in the Lepage-Mackenzie prescription for setting 
the scale. The mass renormalization is small for aM > 3 and becomes large as aM decreases 
below this. Problems in reliably setting the scale result in large uncertainties in determining 
the energy shift and mass renormalization for small aM. These uncertainties underscore the 
need to compute the two-loop corrections to these renormalization parameters. 

This work is partly an extension of previous calculations and is being carried out 
in conjunction with ongoing simulations as part of the NRQCD collaboration |^,^. 

The heavy-quark propagator in lattice NRQCD is briefly described in Sec. 0. The 
energy shift, mass and wave function renormalizations, and the kinetic coupling corrections 



are defined and their calculation is outlined in Sec. |T|. The determination of the Lepage- 
Mackenzie scales is also described in this section. Results of these calculations are presented 
in Sec. |^ A final summary and suggestions for future work are given in Sec. 0. 



II. QUARK PROPAGATOR IN LATTICE NRQCD 

The heavy-quark propag ator g{x) = a^{ij{x)ilj\0)) in lattice NRQCD may be defined 
such that it satisfies an evolution equation given by |^ 

g{x+aei) = Uo^Ulix) [0{x)g{x) + S{x; 0)], 

0(.)^(l-^^y\l-aSH)(l-^-§^y, (1) 
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g{x) = 0, {xi < 0) 

where a is the lattice spacing, n is a positive integer ]![, ip{x) is the heavy-quark field, 
U^{x) is the link variable representing the gauge field along the link between sites x and 
x + aifj, with (e^)^ = 5^^, and uo is a mean-field parameter defined in terms of the average 
plaquette. The nonrelativistic kinetic energy operator Hq is given by 

Ho = , (2) 

and the relativistic and finite-lattice-spacing corrections are included with coupling strengths 
Cj in 

6H = j2c,V„ (3) 

where 

K, = ^(AW.E-E.AW), (6) 

Vi = -^°--(A<±>xE-ExA<±>), (7) 

V. = -4.-B. (8) 

Vt = -gj^{AP),,T. (AWxE-ExA<*>)}. (10) 
1^8 = -^^o- ■ E X E. (11) 

The covariant, tadpole- improved difference operators are defined by 

aA(+) = u^'U, - 1, (12) 
aA(-) = 1 - u,'Ul (13) 
AW = i(AW+A(-)), (14) 



a^'"A(^") = f («oM^. + f/l]-2)'", (15) 
fc=i 



and E and B are the cloverleaf, mean-field-improved chromoelectric and chromomagnetic 
fields, respectively 0. Tree- level removal of the leading cutoff errors in these operators 
produces the operators A'^^\ A'^^^ E, and B. The components of a are the standard Pauli 
spin matrices. The parameter n eliminates doublers and stabilizes the evolution of the quark 
Green's function when n ^ 3/ (Ma). 

The coupling coefficients Cj are functions of g and aM and are determined by matching 
low-energy scattering amplitudes in lattice NRQCD with those from continuum QCD, order- 
by-order in perturbation theory. At tree level, the values of these couplings are unity |0| . The 
0{as) matching determines the radiative corrections to the cj couplings and may possibly 
require the introduction of additional O(a^) -suppressed interactions. 

The operators Vi and V2, although appropriately gauged, will be referred to here as 
kinetic operators because they affect the energy of the free quark in the absence of the gluon 
fields. Their associated couplings ci and C2 will be called kinetic couplings. 

III. KINETIC COUPLINGS AND RENORMALIZATIONS 

The kinetic couplings Ci and C2, the energy shift, and the heavy-quark mass and wave 
function renormalizations may be determined by studying the heavy-quark propagator in 
perturbation theory. An important consideration in the perturbative study of this propa- 
gator is the choice of expansion parameter. Lepage and Mackenzie have advocated a renor- 
malized running QCD coupling defined in terms of the short-distance static-quark potential 
and evaluated at a prescribed mass scale. Their coupling shall be adopted here. 

In this section, Eq. |l]is solved perturbatively to 0{as). The dispersion relation obtained 
from this solution is then matched to that of continuum QCD to O(f^) to determine the 
heavy-quark mass renormalization, the shift in the zero point of energy, and the radiative 
corrections to Ci and 02- The wave function renormalization is then obtained from the 
on-shell residue of the perturbative propagator. Lastly, a brief description of the Lepage- 
Mackenzie coupling av{q*) is presented and their procedure of setting the scale q* is applied 
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to the renormalization parameters and kinetic couplings. 



A. Heavy-quark self-energy 

The heavy-quark self-energy is defined by writing the full inverse quark propagator 
Q^^{p) in the form 

aG-\p)% = Q-\pWS^p - aT.%{p), (16) 

where 2, j are color indices, a, (3 are spin indices, p = (p, P4) is a four-momentum in Euclidean 
space, and Q{p) is the zeroth-order heavy-quark propagator in momentum space given by: 

Q(p) = [e^^^'^-An(p)r„(pf"]"\ (17) 

A (p) = i + ^^!aM!fi + ^)_ ^ (19) 

with Kn(p) = I]j=i sin" (Pi'^/2)- At leading order in perturbation theory, this self-energy is 
given by: 



E(^)(p) + S(^)(p) + S(^)(p) 



(20) 



where T,^^\p) denotes the contribution from the quark-gluon loop diagram shown in 
Fig. [l|(A), T.^^\p) represents the contribution from the tadpole graph shown in Fig. |l](B), 
and T,^'"\p) denotes the sum of the order as contributions from the ci, C2, and link- variable 
renormalization counterterms which arise when one writes Uq = 1 + cXsU^q^ + 0{af) and 
Cj = l + a^cf + 0{al). 

The self-energy S(p) is invariant under interchange of any two spatial momentum com- 
ponents Pi pj and under spatial momentum refiections pj — >■ —Pj, and transforms into its 
complex conjugate under p^ —pi- From these properties, the small-p representation of 
the self-energy may be written: 
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+ hH (EpK) + ...}, (21) 

i=i 

where w = —ip^a. When w is small, the functions fm{w) are represented well by their power 
series expansions fm{w) = Y.'iZo^m 

The moments are obtained from integrals of the form 

where q is the momentum of the gluon in the loop diagrams shown in Fig. |I| and (q) is some 
combination of vertex factors and Feynman propagators. To evaluate these integrals [ITU 



first, a small gluon mass A is introduced to provide an infrared cutoff. Due to singularities 
in the integrands from the zeroth-order quark propagator, it is best to first evaluate the 
integrals over the temporal component of q analytically. Using z = exp(— ig4a), the q^ 
integrations are transformed into contour integrals along the |z| = 1 unit circle which can 
be evaluated by the residue theory. The remaining integrals over the spatial components 
of q are then approximated numerically by discrete sums using aqj = 27mj/N, where the 
rij take all integer values satisfying —N/2 < rij < N/2 for integral A^. The error made 
in this approximation diminishes exponentially fast as ^ cxd, but the rate of decay is 
proportional to A. Convergence in can be dramatically improved by making the following 
change in the variables: qj qj — a sing^, where a = sechw and u satisfies aX = u — tanhw. 
Extrapolation of the results to zero gluon mass is accomplished using a least-squares fit 
to the form: Em=o ^mH'^-^)"' + Em=i ^mH'^'^)"' '^^-^^^ Calculations are performed in the 
Feynman gauge. 

The contributions from to the fm{w) functions appearing in Eq. |2l| can be exactly 

determined. One finds 

(23) 

.(0)/.,A_2.,f2)A 9 9 



/rW = 34'\l + ^ + ^). (24) 



3 ^ V 8MaV n^J AM^a^ 



(2) / ON (2) 

Contributions to the moments fi^-* from the Ci, C2, and mean-field improvement counterterms 
in are then easily determined from the power series expansions of these functions. 



B. Matching 

From the heavy-quark propagator computed to 0{as) in perturbation theory, one finds 
that the on-mass-shell quark satisfies a dispersion relation given to order f ^ by 



r 



where ci;o(p) is the value of w = —ip^a at which the inverse propagator vanishes, Mr is the 
renormalized mass given by = Z^M with = 1 + asi^Q*^ + f^o^'* + ^i^"*), and 

a6uj{p) = Wo + W^ ^5!^ + W2 (Ep>^), (28) 



8M2c. 



Wo = n^^\ (29) 



1^2 = ^^^. (31) 

Note that the rest mass Mr has been removed from the heavy-quark energy in Eq. 

The order as corrections to the heavy-quark propagator in QCD only renormalize the 
quark field and mass to order v'^ so that the dispersion relation for the quark in continuum 
QCD is given by 

p2 (p2)2 

'^QCd(p)-^-^ + ... . (32) 
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If lattice NRQCD is to reproduce the low-energy physical predictions of full QCD, then 

5a; (p) = is required in Eq. ^ that is, Wq = Wi = W2 = 0. Alternatively, since the 

constant term Wq represents an overall energy shift, one could require only Wi = W2 = 

and then simply shift the energies obtained in simulations by an amount asWo/a + Mr 

(2) (2) 

for each heavy quark. Setting Wi = W2 = determines the coefficients and C2 , and 
defining p4 = p4 + iagWo/a, the inverse propagator for small v may then be written: 

ag-\p) ^ Z;^^ (^ip^a + + . . ^ , (33) 

where Z^ = 1 — as , asv"^) is the wave function renormalization. 

A more convenient set of renormalization parameters may be obtained by defining 

p4 = P4 — iasA/a, (34) 
Zm = l + a^B, (35) 
= l + a,(C + Cdiv), (36) 

where C^w = -2{\n a'^ X'^) / (Sir) is the infrared-divergent portion of the wave function renor- 
malization. The parameters A, B, and C can then be calculated using 

B = + ni^^ + nf\ and c = -{nP + n'^p + Cdw). 



C. Expansion parameter av{q*) 

In order to obtain numerical values for the renormalization parameters and the kinetic 
couplings, one must determine the value of the renormalized QCD coupling a^. This can 
be done only after a definition of the running coupling and a procedure for determining the 
relevant mass scale are specified. 

A renormalization scheme, due to Lepage and Mackenzie P], which defines the coupling 
such that the short- distance static-quark potential has no or higher-order corrections is 
particularly attractive. By absorbing higher-order contributions to the static-quark potential 
into the coupling, it is hoped that higher-order contributions to other physical quantities 
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in terms of this coupling will then be small. In this scheme, the running coupling av{q) is 
defined by 

Viq) ^ -^aviq) (37) 

for large q, where Cp = 4/3 is the quark's color Casimir and V{q) is the static-quark 
potential at momentum transfer q. For sufficiently large q, this coupling runs according to 
the usual two-loop relation 



av{q) 



In ln(g''/A 



2 / A 2 N 

v) 



(38) 



where Pq = (11 - |A^ir)/(47r), pi = (102 - fNF)/{16n'^), Ay is the QCD scale parameter, 
and Np is the number of dynamical quark fiavors at mass scale q. 

Lepage and Mackenzie have also devised a simple procedure for determining the scale q* 
at which to evaluate the coupling ay(g*) for a given one-loop process p|. For a one-loop 
perturbative contribution of the form 

h = «v(g*) / ^ m, (39) 

where q is the momentum of the exchanged gluon, they suggest 

Clearly, difficulties with this definition will arise when / d'^q C,{q) ~ 0. Also note that the 
mean value theorem guarantees that q* will satisfy < aq* < 2tt only if ^(g) > for all q 
throughout the region of integration — vr < aq^ < it (or ^(g) < for all such q). 

To evaluate the scales g^, g*(2), and g*(2) corresponding to the gauge-invariant quan- 

(2) (2) 

titles A, B, c\ , and q , respectively, integrals of the form 

= / (0^ In(g^) (41) 



must be calculated (recall Eq. |2^). First, a small gluon mass A is introduced to circumvent 
infrared difficulties. Next, the In(a^g^) factor is approximated by ln{(^(2;) -|- a^(q^ -|- A^)}, 
where 
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C(z) = -15 , (42) 

^ 23^4 + 688^3 + 2358^2 + 688;z + 23 ^ ^ 

and z = exp(— ig4a). Eq. ^ is simply a [2/2] Fade approximant of a'^qf in terms of cos(ag4). 
This introduces a small error in the determination of the scale, but enables an analytical 
treatment of the integration over the temporal component of q. Singularities in the inte- 
grands make such a treatment crucial to reliably and efficiently evaluating these integrals. 
The remaining integrations over the spatial components of q are then approximated numeri- 
cally by discrete sums, as described previously in Sec. ^11 A| , and the A — limits are finally 
taken. 



IV. RESULTS AND DISCUSSION 

Results for A, B, C, cf'\ c^^ and the scales q\, q*^, q*{2), q*{2) are presented in this section. 
Since the wave function renormalization is not a gauge-invariant quantity, the scale corre- 
sponding to C will not be considered here. Choosing a value for aAy from recent simulation 
measurements of ay, numerical estimates of the energy shift, mass renormalization, and the 
radiative corrections to the kinetic couplings are also obtained. 



A. Perturbative coefficients 

(2) (2) 

Values for A, B, C, c\ , and C2 have been determined for a large range of bare quark 
masses and for various values of the stability parameter n. Results were obtained using the 
standard Wilson action Sq^'' for the gluons and a tree-level 0(a2)-improved gluonic action 
Sq^ 10. In addition, two versions of the NRQCD interaction operator 6H were used: one, 
denoted by 6H^^\ containing the full set of interactions Vi — Vg in 5H (see Eqs. P|— P^D; 
the other, denoted by 5H^'^\ including only the 0{v'^) corrections Vi — V5 with no tree- 
level removal of the cutoff errors in V4 and V5 (Eqs. and |^ with the tildes removed). 
Contributions to these parameters from the tadpole-improvement counterterm are given by 
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SB 



6A 




) 



(43) 



(44) 



6C 




(2) 3 



(45) 




18n 



(46) 



Ma + 2n \n 




9 




(47) 



The mean-field correction is Mq = — evr/S, where e = 1 for the simple gluonic action and 
e = 0.78463 for the improved gluonic action. 

The shift A in the zero point of energy, the mass renormalization B, and the kinetic 
coupling corrections Ci and C2 before tadpole improvement (mq = 0) are shown in Figs, y 
and for the simple Wilson gluonic action S^^'^ and the O(f^) version SH^"^^ of NRQCD 
| ri| . Their values after tadpole improvement are shown in Figs. ^ and ^. These quanti- 
ties are gauge invariant and infrared finite. In order to demonstrate the sensitivity of the 
tadpole-improved quantities to the form of the gluon action and the higher-order NRQCD 
interactions, results for B after tadpole improvement using both versions SH^"^^ and 5H'^^^ 
of NRQCD and the simple Wilson 5*^^ and the 0(a^)-improved Sq^ gluonic actions are 
shown in Fig. |^. The wave function renormalization parameter C before and after mean-field 
improvement for Sq^'' and dH^"^^ is presented in Fig. |^. 

Before tadpole improvement, these parameters are all very large, especially as aM de- 
creases. Tadpole contributions dominate and contain power-law divergences which grow as 
aM becomes small. Since high-momentum modes are more strongly damped in the im- 
proved gluon propagator, the ultraviolet-divergent tadpole terms, and hence the couplings 
and renormalization parameters before tadpole improvement, are smaller in the case of the 
improved gluon action. 

Mean-field improvement greatly reduces the magnitudes of all of these coefficients. For 
aM > 5, A is almost independent of aM after tadpole improvement, with values near unity; 
as aM decreases below 5, A also decreases, eventually changing in sign. After tadpole 
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improvement, B has a small negative value at large aM; as aM decreases, B increases, 
becoming positive, reaches a maximum near aM ~ 1, then rapidly decreases, eventually 
falling below zero again. The kinetic coupling coefficients c^i^ and c^2^ vary slowly for aM ^ 
3/2, then rise sharply as aM decreases below this. C also varies slowly, remaining a small 
correction, for aM ^ 3/2 after mean-field improvement. As aM decreases below unity, 
however, C begins to grow quickly in magnitude. For large aM, these tadpole-improved 
parameters are only slightly sensitive to the value of n, the spin- dependent interactions, the 
cutoff improvement of the gluon action, and the 0(f ^) relativity corrections of NRQCD (see 
Fig. ^). As aM decreases, however, sensitivity to these factors grows. The kinetic couplings 
are generally more sensitive to these effects than are the renormalization parameters A, B, 
and C. The mass renormalization parameter B is also very sensitive to the leading kinetic 
interactions Vi and V2. 



B. Scale determinations 

The scales before tadpole improvement (mq = 0) are shown in Figs. M and p| for the simple 



Wilson gluonic action Sq^^ and the O(f^) version 

of NRQCD [ini]. Their values aft 



er 



mean-field improvement are shown in Figs. |l^ and |TT]. The qualitative behavior of these 
scales are not affected by the cutoff improvement of the gluon action and the O(f^) relativity 
corrections of NRQCD. For large aM, the tadpole-improved scales are especially insensitive 
to these factors. As aM decreases, however, quantitative sensitivity grows, becoming rather 
pronounced for q*^2) and q*^2) ■ 

Before tadpole improvement, the scales q\ and g^j are slowly varying functions of aM for 
aM < 5, ranging between 2/a and 3/a in value. As aM becomes large, falls slowly to a 
value just below 2/a, while rises to a value above 5/a. Since the unimproved parameters 
are dominated by highly-ultraviolet tadpole contributions, scales in the range 2/a to 3/a 
are not surprising. The greater damping of the high-momentum modes in the case of the 
improved gluon action Sq^ results in shghtly lower scales in comparison to those obtained 
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using the Wilson action. 

The scales g*(2) and g*(2) also vary slowly between 2/a and 3/a, except in the regions where 

(2) (2) 

the perturbative coefficients and C2 are nearly zero. As these zeros are approached from 
below, the scales tend asymptotically to infinity; approached from above, the scales tend 
to zero. This dramatic behavior is an inevitable consequence of using Eq. ^ to define the 
scale and does not stem from any physical effect; whenever any perturbative coefficient nears 
zero, similar pathological behavior will be observed in the corresponding scale. Such points 
at which the scales tend to zero on one side and asymptotically to infinity on the other shall 
be referred to here as defects. 

For aM > 3, the scale q\ is almost independent of aM with values near 0.8/ a after 
mean-field improvement, revealing a significant reduction from the tadpole removal. As aM 
decreases below 3, q\ falls, becoming small as the zero in A near aM ~ 0.6 is approached. For 
aM below this zero, q\ assumes large defective values which cannot be shown in Fig. |T0|. The 
mass renormalization scale becomes small as the zero in B near aM ~ 0.5 is approached 
from above and the zero near aM ~ 5 is approached from below; it becomes infinitely large 
as one approaches these zeros from the opposite directions. Between these defects, q^ attains 
a maximum value of only about 0.6/a, and as aM becomes large, q^ tends to a value near 
2/a; both of these facts indicate a substantial lowering of q% from tadpole improvement. 

For aM > 3, the scale g*2) is significantly smaller after mean-field improvement and 
varies little; as aM decreases below 3, q*^^) grows quickly, eventually reaching values above 
6/a, indicating an enhancement from tadpole removal. Note that mean-field improvement 
has eliminated all pathological values. q*^2) displays behavior similar to that of g*(2) , except 

Cj C2 

that g*2) grows to much greater values as aM decreases, exceeding 15 /a near aM ~ 1 before 
leveling off. These high values are a result of large, nearly canceling contributions in the 
loop integrals in Eq. ^ 
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C. Size of the corrections 

In addition to the perturbative coefficients A, B, c\ , C2 , and their associated scales, 
the QCD scale parameter aAy must be known before numerical estimates of the first-order 
energy shift aEo = A <yv{(lA)y ^^e mass renormalization Zm — 1 = B and the 

radiative corrections Sci = Ci av{Q*(2)) and 6c2 = C2 ay(5'*{2)) can be determined. The 
most straightforward procedure for measuring aAy is to extract av'(yu) at some scale fi 
from simulations. One possibility is to fit simulation measurements of the logarithm of the 
mean plaquette to its perturbative expansion. Doing this, Lepage and Mackenzie have 
extracted av(3.40/a) from Monte Carlo data for quenched QCD at several values of (3, where 
f3 = Q/g^. At (3 = Q, they find ay(3.40/a) = 0.152, corresponding to aAy = 0.169 in Eq. |3|; 
at /3 = 6.4, they obtain ay(3.40/a) = 0.130, yielding aAy = 0.097. 

The leading-order radiative corrections 6ci and 6c2 to the kinetic couplings are shown 
in Fig. |12| for aAy = 0.169 and using 6H^'^'> and Sq^\ For aM > 1, 6ci and 6c2 are nearly 
independent of aM and adjust ci and C2 above their tree-level strengths by about 10%, a 
small amount. For aM < 1, these corrections grow rapidly as aM decreases, indicating the 
insufficiency of tree-level values for ci and C2 as aM becomes small. These large corrections 
also suggest that the NRQCD approach should be applied cautiously as aM falls below 
unity. 

The defects in q\ and make determination of aEo and Zm — 1 problematical. One 
cannot simply use aEo = A ay{q\) and Zm — 1 = B ay^q]^) near the defects since ay{q) 
is not reliably known for q below and near Ay. Imposing minimum and maximum scales 
also does not work since this produces large, unphysical discontinuities in aEo and Zm — 1 
on crossing the defects. The q* prescription of Eq. ^ fails near the defects and it would 
seem that one must resort to guessing the appropriate scale. One suggestion for doing this 
is to inspect the scale plots and somehow smooth out the scales across the defects. However 
chosen, denote this guessed scale by q to distinguish it from the Lepage-Mackenzie scale 
q*. Estimates of aE^ and Zm — 1 can then be obtained by replacing ay{q*) by ay{q). A 
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better procedure, however, is to replace av{q*) by its expansion in terms of the couphng 
renormahzed at q, namely, 

av{q*; q) = av{q) [l - /So In {q*'^/f) ay (g)] . (48) 

This incorporates information contained in q*, yet yields results which are continuous across 
the defects. 

A range of estimates for the energy shift aEo and mass renormalization Zm — 1 using 



Eq. ^ are presented in Fig. |T3|. Energy shift estimates A av{qA', qA) are shown choosing = 
0.8/a and 0.6/a for all aM values and using aAy = 0.169 in a;y(g^). Mass renormalization 
estimates B Oiv{q*B]qB) are shown for the choices = 1.0/a and 0.6/a for all aM. These 
choices for qA and qs are based on an examination of the scales in Figs. ^ and [ij. The 
energy shift aE^ is remarkably independent of aM, in contrast to the perturbative coefficient 
A which falls to below zero as aM decreases. Sensitivity to the choice of qA grows as aM 
decreases. The mass renormalization is very small for aM > 3 and is only mildly sensitive 
to the choice of qs', — 1 becomes zero near where the perturbative coefficient B vanishes. 
As aM decreases below this, however, — 1 grows quickly and continues to grow even after 
B has fallen to below zero. As aM becomes small, estimates of the mass renormalization 
become very large and the sensitivity of these estimates to g^ becomes pronounced. Both 
of these facts and the sensitivity of the energy shift to (Ja for small aM underscore the need 
to compute the two-loop corrections to Z^ — 1 and aEQ. 



V. CONCLUSION 



The heavy-quark propagator in NRQCD was calculated to 0{as) in tadpole-improved 
perturbation theory. The mass renormalization, shift in the zero point of energy, and the 
radiative corrections to the two kinetic couplings Ci and C2 were determined by matching 
the on-mass-shell dispersion relation for the heavy quark obtained from the NRQCD quark 
propagator to that of continuum QCD to O^v"^). In addition, the scales at which to evaluate 
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the running coupling for these quantities were obtained for the first time using the Lepage- 
Mackenzie prescription. Defects in these scales were found. The wave function renormaliza- 
tion was extracted from the on-shell residue of the perturbative propagator. Results were 
obtained using two different gluon actions Sq^^ and Sq^ and two versions of NRQCD, 5H^'^^ 
and 5if , differing only in the order of the relativity and finite- lattice-spacing corrections 
retained. 

Using a typical value for oAqcd from recent simulation measurements of ag, numerical 
estimates of the energy shift, mass renormalization, and the radiative corrections to ci and 
C2 were also obtained. The radiative corrections to the kinetic couplings were found to be 
small for aM > 1, being roughly 10% of the tree- level couphng strengths, and to grow 
quickly as aM decreased below unity. Problems in rehably setting the scale complicated 
the determinations of the energy shift and mass renormalization for some values of aM. A 
simple modification of the running coupling was necessary to obtain a range of estimates for 
these parameters near these aM values. The mass renormalization was found to be small 
for aM > 3 and to become large as aM decreased. Large uncertainties in determining the 
energy shift and mass renormalization for small aM demonstrated the need to compute the 
two-loop corrections to these quantities. 

Future work includes the evaluation of the other NRQCD couplings and the two-loop con- 
tributions to the mass renormalization and the energy shift. Computation of the couplings 
C3 — C7 has already begun. 
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FIGURES 




FIG. 1. Two Feynman diagrams which contribute to the heavy-quark self-energy. A curly line 
denotes a gluon; a double solid line denotes a heavy quark. 
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FIG. 2. The energy shift parameter A, heavy-quark mass renormahzation parameter B, and 
the kinetic coupUng coefficients c^^ and c^* before tadpole improvement (uq = 1) against the 
product of the bare heavy-quark mass M and the lattice spacing a using S^P and 5if(2)^ The 
dot-dashed curve is the solid curve shows 5, the dashed curve is (W and the dotted curve 
indicates C2 . For aM > 3, the stability parameter n is set to unity; for 1.5 < aM < 3, n = 2 is 
used; for 1 < aM < 1.5, n = 3 is used; and for 0.5 < aM < 1, n = 6 is used. 
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FIG. 3. The energy shift parameter A, heavy-quark mass renormahzation parameter B, and 
the kinetic couphng coefficients c^^ and before tadpole improvement (uq = 1) against 1/aM 
usmg S'Jf^ and 5H(^1 The dot-dashed curve is A, the sohd curve shows B, the dashed curve is 
c^^'* , and the dotted curve indicates c^^ . The stabihty parameter n is set to unity. 
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FIG. 4. The shift A in the zero point of energy, the heavy-quark mass renormahzation parame- 
ter i?, and the kinetic couphng coefficients c\ and C2 a/ter tadpole improvement (no = 1 — as7r/3) 
against the product of the bare heavy-quark mass M and the lattice spacing a using and 
5/^(2)^ The dot-dashed curve is the solid curve is the dashed curve is Sy \ and the dotted 
curve indicates S^^ ■ The stability parameter n assumes the same values as used in Fig. ^. 
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FIG. 5. The shift A in the zero point of energy, the heavy-quark mass renormahzation parame- 
ter B, and the kinetic couphng coefficients c\ and C2 a/ter tadpole improvement {uq = 1 — a^Tr/S) 
against 1/aM using S^^^ and The dot-dashed curve is A, the solid curve is B, the dashed 



curve is c^^ , and the dotted curve indicates C2' ■ The stability parameter n is set to unity. 
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FIG. 6. The heavy-quark mass renormahzation parameter B against the product of the bare 
heavy-quark mass M and the lattice spacing a. Results using the improved gluon action S'^'* and 
the full set of NRQCD interactions SH^"^^ are shown solid curve. Results using and 5/7(2) 
are shown as a dot-dashed curve. The dashed curve indicates the results obtained using the simple 
gluon action Sq^^ with 5H^^\ while the dotted curve shows the results using 5jf ) with 5H^^). 
The stability parameter n assumes the same values as used in Fig. ^ 
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FIG. 7. The heavy-quark wave function renormalization parameter C against the product of the 
bare heavy-quark mass M and the lattice spacing a using ^g^^ and 5^(2). The dashed curve shows 
C before tadpole improvement (uq = 1) and the solid curve shows C after mean-field improvement 
(uq = 1 — as7r/3). The stability parameter n assumes the same values as used in Fig. ^. 
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FIG. 8. The energy shift scale heavy-quark mass renormahzation scale glj, and the kinetic 
coupling coefficient scales q*,2) and q*,2) before tadpole improvement (uq = 1) against the product 
of the bare heavy-quark mass M and the lattice spacing a using S^^^ and (5/7(2). dot-dashed 
curve is ag^, the solid curve shows ag^, the dashed curve is ag*(2), and the dotted curve indicates 
ag*(2) • The stability parameter n assumes the same values as used in Fig. ^. 
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FIG. 9. The energy shift scale q\, heavy-quark mass renormalization scale g^, and the kinetic 
coupling coefficient scales q*(^2) and q*(^2) before tadpole improvement {uq = 1) against 1/aM using 

Cj C2 

S'^^^ and (5i?(2)_ The dot-dashed curve is aq\, the solid curve shows aq*^, the dashed curve is aq*^^), 

^1 

and the dotted curve indicates ag*(2) . The stability parameter n is unity. 
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FIG. 10. The energy shift scale q\, heavy-quark mass renormahzation scale g^, and the kinetic 
coupling coefficient scales q*^2) and q*^2) a/ier tadpole improvement (uq = 1 — a^vr/S) against the 
product of the bare heavy-quark mass M and the lattice spacing a using and <5F(2). The 

dot-dashed curve is aq\, the solid curve shows ag^, the dashed curve is ag*(2), and the dotted curve 
indicates aq*^^2) ■ The stability parameter n assumes the same values as used in Fig. |2[ 
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FIG. 11. The energy shift scale q\, heavy-quark mass renormalization scale q^, and the kinetic 
coupling coefficient scales q*(^2) and g*(2) a/ter tadpole improvement {uq = 1 — a^Tr/S) against 1/aM 
usmg S^J^^ and The dot-dashed curve is aq^, the solid curve shows aq^, the dashed curve 

is ag*(2) , and the dotted curve indicates og* 2) . The stability parameter n is unity. 
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FIG. 12. The first-order radiative corrections to the kinetic couplings ci and C2 against aM 
using and 6H^^\ The dashed curve is 6ci = c^^ ay(9*(2)) and the dotted curve indicates 

6c2 = C2 ctv{Q*(2))j using aAv = 0.169. The stability parameter n assumes the same values as 
used in Fig. 0. 
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FIG. 13. Estimates of the energy shift uEq and mass renormahzation — 1 against aM using 
Sq^^ and 5H^'^\ The dot-dashed curves labeled a and b are energy shift estimates A av{q\;qA) 
for qA = 0.8/a and 0.6/a, respectively, using aKy = 0.169 in av^qA)', the unlabeled dot-dashed 
curve is A ay{q\). The solid curves labeled c and d show B ayiq^i qs) for qB = l-O/o and 0.6/a, 
respectively; the unlabeled solid curve is B av{q*B)- Values for the stability parameter n are the 
same as in Fig. 0. 
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